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I We study the quantum-mechanical tunneling phenomenon in models 

Q, 

which include the existence of non-equivalent vacua. For such a purpose 
we evaluate the euclidean propagator between two minima of the potential at 



! issue in terms of the quadratic fluctuations over the corresponding instantons. 



o 

(N 
O 



The effect of the multi-instanton configurations are included by means of the 
alternate dilute-gas approximation. 



Ph! I- INTRODUCTION, 

g 

The tunneling through classically forbidden regions represents one of the most striking 



phenomenon in quantum theory and therefore plays a central role in many areas of modern 



^ ■ physics. On the other hand, together with the operator formalism of quantum mechanics we 



have an equivalent description by means of path-integrals. In such a case the Schrodinger's 
equation is substituted by a global approach where the quantum mechanical time evolution 
is analysed in terms of a functional integration. Qualitatively speaking, the path-integral 
representation corresponds to a sum over all histories allowed to the physical system we 
are dealing with. To be precise, we need to take into account an imaginary exponential of 
the classical action associated with every path which fulfill the appropiate initial and end 
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points conditions. Of course, the quantum amplitude so built is difficult to handle due to 
the oscillating character of the exponential at issue. To avoid problems of this sort we carry 
out the change t — > —it (known in the literature as the Wick rotation). In doing so, we can 
take advantage of the euclidean version of the path-integral which represents by itself a new 
tool for describing relevant aspects of the quantum theory. 

Almost from the very beginning of the subject a semiclassical treatment of the tunneling 
phenomena (ranging from periodic-potentials in quantum mechanics to Yang-Mills models in 
field theories) has been performed by means of the so-called instantons 0. Going to more 
physical terms, the instantons represent localised finite-action solutions of the euclidean 
equation of motion. To be specific, the euclidean equation is the same as the usual one for our 
particle in real time except that the potential is turned upside down. Although more massive 
than the perturbative excitations, the instantons themselves become stable since an infinite 
barrier separates them from the ordinary sector of the model. The stability is reinforced by 
the existence of a topological conserved charge which does not arise by Noether's theorem 
in terms of a well-behaved symmetry, but characterizes the global behaviour of the system 
when the imaginary time is large enough. Accordingly, it comes by no surprise that these 
classical solutions have been considered in the literature under the name of topological 
configurations. Once the appropiate classical solution is well-known, we make an expansion 
around the topological background to evaluate the quadratic fluctuations which arise in 
terms of the functional determinant of a second order differential operator. The integration 
is solved within the gaussian scheme except for the zero-modes which appear due to the 
invariances of the model. To deal with these excitations we introduce collective coordinates 
so that ultimately the gaussian integration is carried out only along the directions orthogonal 
to the zero-modes. As a functional determinant includes an infinite product of eigenvalues, 
the result should be in principle a highly divergent expression. Fortunately we can regularize 
the fluctuation factors by means of the ratio of determinants. 

To start from scratch let us describe the instanton calculus for the one-dimensional 
particle as can be found in f2|. The reader finds there a comprehensive description of this 
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subject. To be specific we assume that our particle moves under the action of a confining 
potential V(x) which yields a pure discrete spectrum of energy eigenvalues. If the particle 
is located at the initial time U = —T/2 at the point X{ while one finds it when tf = T/2 
at the point Xf, the functional version of the non-relativistic quantum mechanics allows 
us to express the transition amplitude in terms of a sum over all paths joining the world 
points with coordinates (—T/2,Xi) and (T/2,Xf). Making the change t — > —ir, known in 
the literature as the Wick rotation, the euclidean formulation of the path-integral reads 

< x f \ exp(-HT)\ Xl >= N(T) J [dx] exp{-S e [x(T)}} (1) 

where H represents the hamiltonian of the model, the factor N(T) serves to normalize the 
amplitude while the symbol [dx] indicates the integration over all functions which fulfil the 
boundary conditions at issue. Now we resort to the euclidean action S e , i.e. 
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dr (2) 



r T/2 
-T/2 



/ 

J-l 



whenever the mass of the particle is set equal to unity for notational simplicity. In this 
approximation we take for granted the existence of a x c (t) configuration which represents a 
stationary point of the euclidean action. Next we perform the expansion of a generic x(t) 
according to 

x(r) = x c (r) +J2cj xjir) (3) 
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where as usual Xj(r) stand for a complete set of orthonormal functions 

rT/2 

J Xj{r) x k (r) dr = S jk (4) 



vanishing at our boundary, i.e. Xj(±T/2) = 0. The eigenf unctions Xj(r) (with eigenvalues 
€j) appear associated with the so-called stability equation given by 

-^M + V tt [x e (r)]x j (r) = e j x j (r) (5) 
To sum up, the euclidean transition amplitude reduces itself to 



< x f \exp(-HT)\xi >= N(T) exp(—S eo ) JJ^ 



1/2 



(6) 



where S eo represents the classical action associated with the configuration x c (t) while the 
product of eigenvalues is usually written as 



-1/2 



Det 



-1/2 



(7) 



according to the conventional notation of the finite-dimensional case. To fix the factor N(T) 
we resort to a well-known problem where 



v 



v( x ) = y x 



(8) 



so that V"(x = 0) = v 2 . As corresponds to the harmonic oscillator the relevant amplitude 
is given by 



< x f = 0| exp(-H ho T)\xi = >= N(T) ^Det 



d 2 



dr 2 



-1/2 



(9) 



Next the evaluation of (^) is possible according to the general method exposed in . In 
other words 



<x f = 0\ exp(-H ho T)\xi = 0>- 



v 



71 



1/2 



(2 sinh uT) 



-1/2 



(10) 



II. THE SEXTIC-MODEL. 



As anticipated in the introduction the tunneling phenomenon appears as first in the 
study of systems which have two or more degenerate classical minima separated by potential 
energy barriers. The standard example would be the point-particle in a double-well potential 
where, because of tunneling precisely, the two first eigenstates are the spatially even and odd 
combinations of harmonic oscillators centered at the bottom of the respective wells. In more 
physical terms, the degeneracy of the energy eigenvalues is broken so that the splitting term 
results proportional to the barrier-penetration factor. The euclidean configuration which 
leads the process is of course the instanton whose stability properties are decided by the 



behaviour of the small fluctuations in its neighborhood. As a matter of fact the stability 
is assured by a clever combination of translational invariance and conserved topological 



From a classical point of view we find three minima located at 1 and 

x Q = 0. Two of them, namely x_ and x + , are equivalent since can be connected by means of 
the discrete symmetry x — > — x at issue. However the third minimum x Q is invariant under 
the action of such a symmetry. In other words, the central vacuum is not identical with the 
others located at both sides. When considering the limit uo 2 ^> 1 the energy barriers are high 
enough to decompose the system into a sum of independent harmonic oscillators. However 
the existence of finite barriers between the different wells of the potential yields a relevant 
tunneling phenomenon so that ultimately the symmetry x — > —x is not spontaneously 
broken at quantum level and the expectation value of the coordinate x computed for the 
ground-state is zero as corresponds to the even character of the potential V(x). 



The question we wish to address at this point should be the explicit description of 
the tunneling in the euclidean version of the path-integral. As regards the one-instanton 
amplitude we take into account the transition amplitude between x Q = and x + = 1. For 
such a purpose we need a classical configuration with Xi = at U = —T/2 while Xf — 1 
when tf = T/2. It is customarily assumed that T — > oo mainly because the explicit solution 
of the problem is much more complicated for finite T. 

First of all we can find the explicit form of the instanton x c (r) just by integration of a 
first-order differential equation, i.e. 



charges. In the following we consider the triple-well potential V(x) given by [|j] 




(11) 



A. The one-instanton contribution. 




(12) 
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where we recognize the quantum mechanical version of the Bogomol'nyi condition [Q. Now 
we solve (0) by a simple quadrature so that 



x c (t) = \J[1 + tanhu;(r — r c )]/2 



(13) 



where as usual the parameter r c indicates the point at which the instanton makes the jump. 
As expected equivalent solutions are obtained by means of the transformation r — > — r and 
%c( T ) ~~ > —x c {t) so that adjoint minima of the potential can be connected by means of a 
topological solution in a systematic way. In addition we have that S eo = to /4. On the 
other hand we need classical configurations for which x- = and x + — 1 at large but 
finite values r = ±T/2. However the explicit form of the instantons that appear in the 
literature corresponds to infinite euclidean time. Fortunately the difference is so small that 
can be ignored mainly because we are interested in the limit T — > oo. Our description of 
the one-instanton amplitude between x_ = and x + = 1 takes over 



<x f = l\ exp(-HT)\x, t = >= N(T) ^Det 



d 2 



+ 1/ 



-1/2 



Det [-{d 2 /dT 2 ) + V"[x c {r)]] 



-1/2 



exp(-S' eo ) 



(14) 



[ Det[-(d 2 /dr 2 ) +u 2 ) J 

where we have multiplied and divided by the determinant associated with the harmonic 
oscillator of frequency v. As regards the determinant built over the instanton itself we find 
a zero-mode x q (t) which could jeopardize the computation procedure as a whole. However 
this eigenvalue e a = comes by no surprise since it reflects the translational invariance of 



the system. We can discover the existence of this zero- mode starting from (13). Including 
the adequate normalization one can check that 

, i 1 dX r 



XAT) 



(15) 



V ' S eo dT 

is just the solution of (|5|) with e G = 0. The way out of this apparent cul-de-sac is simple. 
The integration over c Q (see (§)) becomes equivalent to the integration over the center of the 
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instanton r c . To fix the jacobian of the transformation at issue we consider a first change 
like 



Ax(t) = x (t) Ac (16) 

Going back to the general expression written in (§) we find that under a shift Ar c the 
global effect corresponds to 

Ax(t) = -y/sZ x {t)At c (17) 

Finally the identification between ([16]) and (|T7| ) yields 

dc Q = \fsZ dr c (18) 

where the minus sign disappears since what matters is the modulus of the jacobian at issue. 
To sum up we have that 



j Det [-(d 2 /dT 2 ) + V"[x c (r)}} 
{ Det[-{d 2 /dr 2 ) + u 2 ] 



-1/2 



j Det'[-(d 2 /dr 2 ) + V"[x c (T)]] }- 1/2 fsZ ,. Q , 
\ Det[-{d 2 /dr 2 ) + v 2 ] J V27T c { ' 

where Det' stands for the well known reduced determinant once the zero-mode has been 
removed. To make an explicit computation of the quotient of determinants we resort to the 
Gelfand-Yaglom method where only the knowledge of the large-r behaviour of the classical 
solution x c (t) is necessary |J. If O and P represent a couple of second order differential 
operators, whose eigenfunctions vanish at the boundary, the quotient of determinants is 
given in terms of the respective zero-energy solutions / (r) and g (T~) according to 

DetO f (T/2) 



DetP 9o(T/2) 
whenever the eigenfunctions at issue satisfy the initial conditions 



(20) 



f o (-T/2) = g o (-T/2) = 0, ^(-T/2) = ^(-T/2) = 1 (21) 



For the zero- mode g (r) associated with the generic harmonic oscillator of frequency v 
we have 

9o(t) = ~ sinh[z/(r + T/2)} (22) 

so that now we need the explicit form of the solution f (r) which corresponds to the topo- 
logical configuration written in (|13|). Starting from the x (t) zero-mode we write a second 
solution y (r) given by 

y (r)=x (r) [-^ (23) 

In doing so we may summarize the asymptotic behaviour of x q (t) and y {^) as follows 

{Cexp(— 2ut) if r — > oo 
(24) 
Dexp(uT) if r — > — oo 

exp(2tjr) /4uC if r — > oo 
y (r) ~ <( (25) 
- exp(— ujt)/2ujD if t — > — oo 



where the constants C and D derive from the explicit form of the derivative of (|13"|). Starting 
now from the linear combination of x G (r) and ?/ (r) given by 

/ (r) = Ax (r) + By (r) (26) 

the employ of the initial conditions leads us to 

Ur) = x (-T/2)y (r) - y (-T/2)x (T) (27) 

From this expression we can extract the asymptotic behaviour of / (t), i.e. 

fo(T/2) ~ ^ exp(ujT/2) if T ^ oo (28) 

Now we need to take into account the lowest eigenvalue of the stability equation to 
obtain the value of the ratio of determinants. We can explain the situation as follows: the 
derivative of the topological solution does not quite satisfy the boundary conditions for the 
interval (— T/2, T/2). When enforcing such a behaviour, the eigenstate is compressed and 
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the energy shifted slightly upwards. In such a case the zero-mode x (t) is substituted for 
the /a(t) which stands for 

d 2 fx(r) 



dr 2 

whenever 



+ V"[x c {r)\h{r) = AA(r) (29) 



U-T/2) = f x (T/2) = (30) 



Resorting now to the lowest order in perturbation theory we get 

h(r) ~ fo(r) + A ^ 



(31) 



A=0 



so that 



h(r) = f (r) + A f [x (r)y (s) - y (r)x (s)] f (s) ds (32) 

J-T/2 

The asymptotic behaviour of /°(t), x (r) and y G (r), together with the condition 
fx(T/2) = allows us to find this lowest eigenvalue A, i.e. 

A = 2ujD 2 exp(-cuT) (33) 

The Gelfand-Yaglom method provides us with the final expression of the quotient of 
determinants whenever we choose for the frequency v of the harmonic oscillator of reference 
the average of the frequencies of the central and lateral wells. In other words u = 3u/2. We 
can compare this situation with the well-grounded double-well model where the two minima 
of the potential are equivalent so that the aforementioned average is out of order. With 
this information we can write the one-instanton amplitude between the points X{ = and 
Xf = 1, i.e. 

<x f = l\ exp{-HT)\xi = >= 

o x 1/2 HT 

—J (2sinh3cuT/2)- 1/2 sj S eo J— exp(-S eo ) to dr c (34) 



As expected we get a transition amplitude just depending on the point r c at which the 
instanton makes precisely the jump. Next we should take into account the configurations 
constructed out of instantons and antiinstantons which mimic the behaviour of a trajectory 
strictly derived from the euclidean equation of motion. In doing so we get an additional 
bonus since the integration over the centers of the string of instantons and antiinstantons is 
carried out in a systematic way. But this point is the subject of the next subsection. 



B. The dilute-gas approximation. 

Since the above calculations were carried out over a single instanton, it remains to identify 
the contributions associated with a string of widely separated instantons and anti-instantons 
along the r axis. It is customarily assumed that these combinations of topological solutions 
represent no strong deviations of the trajectories just derived from the euclidean equation of 
motion without any kind of approximation. We evaluate the functional integral by summing 
over all such configurations, with k instantons and anti-instantons centered at points n, Tfc 
whenever 

T T 
~- < 7-1 < ... < Tfc < - (35) 

If the regions where the instantons (anti-instantons) make the jump are narrow enough, 
the action of the proposed path is almost extremal. In this approach the total action is 
given by the sum of the k individual actions. This scheme is well-known in the literature 
where it appears with the name of dilute-gas approximation [0. The translational degrees 
of freedom of the separated k topological configurations yield an integral of the form 

fT/2 rr k rr 2 (u)T) k 

/ ujdr k / udr k -i... / codri = — -— (36) 

J-T/2 J-T/2 J-T/2 k\ 

As regards the quadratic fluctuations we have now that the single ratio of determinants 
transforms into 

f^V /2 ,0^0.^/^-1/2 \Det>[-{d*/dT*) + V"[x c {T)]}\- 1/2 



( — (2 sinh 3cjT/2V , , 
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2tt/ 



(37) 



\ Det[-(d 2 /dr 2 ) + 9u 2 /4] J 
according to the limit of the factor associated with the harmonic oscillator when T is large. 
When going to the dilute-gas approximation one considers a set of instantons and anti- 
instantons so that each topological configuration starts where its predecessor ends. For 
our amplitude the total number k of instantons plus anti-instantons must be odd. The 
combinatorial factor F associated with the number of possible configurations corresponds 
to F = 2 <yk ~ 1 ^ 2 because the closed paths starting and coming back to the point x Q = 
require in a systematic way instanton-anti-instanton pairs. We remind the difference with 
the double-well potential where the instantons strictly alternate with the anti-instantons 
since the problem has only two minima and the combinatorial factor is not necessary. Now 
we can express the complete transition amplitudes for the triple-well potential so that 

<x f = l\eM-HT)\ Xi = >= (-) exp(-3o;T/4)X: ^-^ (38) 

where as usual d stands for the instanton density, i.e. 




d = \hr~ V s ™ exp(-S eo ) (39) 

37T 



To sum up 



'3uA 1/2 

< x f = 1| exp(-HT)\x t = >= ( — exp(-3cjT/4) sinh(ujTd) (40) 



Ait J 

If we return to the euclidean transition amplitude written in ([!]), the insertion of the 
pure discrete spectrum of energy eigenfunctions, namely 

H\n >= E n \n > (41) 

allows us to write that 

<Xf = l\ exp(-HT)\xi = >=J2 ex P(~ E nT) <x f = l\n>< n\x t = > (42) 

n 

Next we denote by E Q , E\ and E 2 the energies of the ground-state and the two first 
excited levels of our problem. As the triple-well potential we are dealing with is even we 
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know that < l\xi = > vanishes so that the limit T — > oo in ( |40| ) provides us with the 
following energy eigenvalues 

^ 3u; , _ 3a; „ 3co> . 
E = — -tud, E x = —, E 2 = — + ud (43) 

The tunneling transfers along the whole real axis the gaussian wave functions constructed 
over the points x_ = — 1, x + = 1 and x Q = so that the splitting term results proportional 
to the barrier-penetration factor. In addition the average of the harmonic frequencies over 
the non-equivalent minima of the potential serves as the central position for the splitting. 



REFERENCES 

[1] A. Polyakov, Nucl. Phys. B120 (1977) 429. 

[2] H. Kleinert, Paths Integrals in Quantum Mechanics, Statistics and Polymer Physics. Singapore: 
World Scientific (1990). 

[3] P. Ramond, Field Theory: A Modern Primer. New York: Addison- Wesley Publishing Company 
(1992). 

[4] J. Casahorran, Phys. Lett. A283 (2001) 285. 

[5] E. B. Bogomol'nyi, Sov. Journ. Phys. 24 (1973) 1888. 

[6] I.M. Gelfand and A.M. Yaglom, Journ. Math. Phys. 1 (1960) 48. 

[7] S. Coleman, Uses of Instantons in The Whys of Subnuclear Physics. Ed. A. Zichichi. New York: 
Plenum Press (1979). 



12 



